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LI.  LNTRODUCTION 


Ae OVERVIEW 


Sputtering is the removal of surface atoms due to energetic 
particle bombardment. It is caused by collisions between incoming 
particles and atoms in the selvage, iee., the near surface of solid 
material. First discovered by Grove [1] in 1853 and Faraday in 1854, 
it took about 50 years until the physical process involved ws 
recognized and another 50 years until a quantitative theory began to 
develop. Grove [1] noticed the disintegration of cathodes in glow- 
discharge tubes and observed, that the cathode material was deposited 
on the surrounding glass walls. He called this process “cathode 
Sputtering”. 50 years later, Goldstein [2] presented evidence that the 
Sputtering effect was caused by positive ions of the discharge hitting 
the cathode. 

Stark [3,4] introduced a so-called ‘hot-spot' model of sputtering 
(now known as ‘evaporation theory') and later a collision theory, which 
viewed sputtering as a sequence of binary collisions initiated by the 
bombarding ion. Stark found that at low energies the yield, Y, 
increases with increasing ion energy, and interpreted this result on 
the assumption that an increasing amount of energy is transferred to 
the target atoms near the surface. At higher energies, increasing 
penetration depth of the ion causes a relatively smaller sputtering 
yield and thus a plateau in the yield curve. The reason for this is 
a decrease in the “size” of the ion as its energy increases. Stark 


considered his hot-spot model, the evaporation of target mterial from 





a microscopically small region with high local temperature, and his 
collision theory as two different views of the same process. 

Von Hippel and Blechschmidt [5], who developed a quantitative form of 
Stark's hot spot model, made the first attempt to formulate a 
Sputtering theory on the basis of local heating. Subsequent investi- 
gators found that this model is inadquate. Wehner [6] showed that 
local evaporation alone couldn't explain sputtering phenomena from 
single crystals. In 1934 Lamar and Compton [7] suggested that binary 
collision processes dominate in light-ion sputtering and local 
evaporation is dominant in heavy-ion sputtering. Their ideas led to 
the 'thermal spike’ concept, in which a long lived high temperature 
volume presists in the target after the collision cascade is complete. 
Wehner's observations revived the interest in the collision theory. 
Keywell [8] applied neutron transport theory to Stark's collision 
model. Due to his and subsequent calculations by Harrison [9], 
probability concepts and collision cross sections were introduced as 


important parts of the sputtering theory. 


Be PHYSICAL MEANING 

Sputtering mechanisms operate in the near surface region of the 
target, and are caused by a collision cascade. The processes involved 
are simular to those causing radiation damage in the bulk mterial of a 
solid. Sputtering happens far from thermal equilibrium, which means 
that it isn't to be confused with evaporation. Generally an incident 
particle collides with atoms of the bulk material and transfers energy. 


If the transffered energy is larger than the cohesive energy at the 





lattice site a primary recoil atom is created- This recoil atom 
collides with other target atoms and distributes the energy into a 
collision cascade. The condition for an atom to be ejected is that the 
kinetic energy associated with the normal component of its velocity is 
larger than the potential energy binding it to the surface. The 
sputtering yield is defined as the average number of target atoms 


ejected per incident particle. 


C. LATER WORK 

Although considerable discrepancies exist in the experimental 
results produced by various scientific groups, standards have developed 
which define the basic requirements for reproducible and reliable 
experiments. This development is shown in reviews of Behrisch [10], 
daaneety [11], Carter and Colligen [12], and Pleshivtsev [13]. Over 
the last 30 years theoretical efforts have been increasingly successful 
in explaining the main features of sputtering. 

The analytic theoretical work in this area was done by Sigmund 
[14]. He presented a systematic theory of sputtering on the basis of a 
minimum of assumptions. Basic assumptions concern the characteristics 
of a single collision (differential cross section, inelastic Stopping), 
the structure of the medium (random or crystalline), the structure and 
binding forces of the surface, and some factors that depend on specific 
experimental arrangements like beam-target geometry, high or low ion 
energy, etce The yield calculation consists of the following Steps: 

1) to determine the amount of energy deposited by energetic 


particles (ion and recoil atoms) near the surface; 





2) to convert this energy into a number of low-energy recoil atoms; 
3) to determine how many of these atoms come to the surface; 
4) to select those atoms that have enough energy to overcome the 
surface binding forces. 
Considering, for example, the backsputtering yield from a plane 
surface, the outward current of the target atoms through the surface 


plane at x = QO is found to be: 





, r dE d°2, 
J(E 2 )dE dQ = YFL(E, @, 0) oe | cos QI 
- os E Ide /dx| Ai 


Where E is the energy of the ion impinging at an angle O onto a plane 
surface, Fp is the deposited energy density, Y is the number of 
projectiles per unit time, 2, is the unit solid angle, and [f, is a 
parameter which depends on ¥ and the scattering cross section. The 
sputtering yield is found by integrating this equation over E, and 
©, and dividing the sputtered current by the incident current ¥. 
This leads to 

Y = A Fgte, @, 9) 


with 


_T 7+ 4 
A= mn ii O if d (cos 0) Icos QO, ! P (Eo> 0.) 


2  €£ IdE /dx] 
O Oo 


where P(E,, 95) is the probability that an atom will escape from 
the surface. The simplest model for the surface binding of a metal is 
based on a planar surface barrier U,- In this case, 


1 for Ey cos70, ym 


men. 50.) 4 
a> QO for E, cos“0, < Us 
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which leads to a simple form: 


ide 
an? NCU 
QO O 
with 
aT 2 
“0 = = “0 * BM? 


Where A, = 24 and apy = 0.219 A. 


De. APPLICATIONS OF SPUTTERING 

For a long time sputtering was regarded as an undesired effect 
which destroyes the cathodes and grids in discharge tubes. In fusion 
research, plasma contamination by metal atoms is still one of the mijor 
problems. Sputtering also may cause the distruction of diaphragms and 
targets in accelerators and high-voltage electron microscopes. 

On the other hand, sputtering has become an indispensible process 
in modern technology. It allows the removal of very tightly bound 
surface layers ona nearly atomic scale. Sputtering is used in 
Sputter-ion sources, for obtaining reproducibly clean surfaces and, for 
example, in micro-machining and depth profiling of thin films. 

Sputtering also can be used to obtain information about 
processes occuring on a atomic scale. Many scientists, including 
Kingdon and Langmuir (15), Winters (16), and Taglauer and Heiland (17), 
Studied the sputtering of chemisorbed gas to investigate mechanisms 
Which lead to the emission of multi-component mterials. This thesis 


uses a computer Simulation to varify one assumption made by Winters. 
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E. COMPUTER SIMULATIONS 

Due to their high operating speed, computers are attractive for 
calculations concerning sputtering events. A computer simulation can 
never replace an experiment, but by matching the results, a simulation 
can help to understand and identify sputtering mechanisms, and to 
interpret experimental results. A computer simulation is not a theory, 
it is a mathematical tool used to test the fundamentals of a theory 
{19}. 

Computer Simulations can be divided in two categories: time-step 
and event-store models... Event~store programs move from event to event 
skipping the intervening time. For sputtering, they are based on 
binary collision theory. These fast, complex programs are useful when 
the physical model is well understood. The program used at the Nval 
Postgraduate School is a time-step model based on simultaneous multiple 
interactions (MI). MI simulations of physical systems usually require 
the calculation of about 100 individual ion trajectories (cascades) , 
while binary collision simulations use several thousands. This, and 
the fact that in ML simulations mechanisms are easy to trace, makes the 
time~step model attractive. 

This thesis is based on QRADS, a refined version of QRAD, the 
standard program in use at the Naval Postgraduate School for some time, 
Which is a full-lattice simulation computer program, designed to model 
the dissipation of an incident ion's momentum ina single crystal 
target, using classical mechanics. The program has the ability to 
generate different surfaces of several crystal structures. To obtain 


reacted Surfaces, a special subroutine places adatoms on the surface. 
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After the lattice is build, the potential functions are calculated. 

The impotant mathematical manipulations, concerning the trajectory, are 
done by the subroutine INTEG (integrator). This subroutine also tests 
for the termination of the trajectory. Finally, subroutine LEGIT 


determines Which atoms Were actually ejected. 


Il. OBJECTIVES 


Some experimental results exist for sputtering from 
polycrystalline BCC structures [19] but only a few for clean single 
crystal surfaces or chemically reacted single crystal surfaces [19]. 
Winters [18] investigated the sputtering cross sections of tungsten and 
molybdenum bombarded by argon and xenon ions in the energy range from 
0.3 keV up to 5.0 keV. Figure 1 shows the cross sections as a function 
of incident ion energy. Winters' experiments show that the sputtering 
yieid of nitrogen tends to increase as the atomic weight of the 
surrounding atoms increases. He furthermore concludes, that the 
difference in the cross sections for the sputtering on nitrogen from 
WCO0O1) and Mo(001) is dominated by the absorbate-substrate mass 
difference since the other important parameters are similar. 

The primary objective of this thesis is to study Sputtering from 


nitrogen reacted (001) tungsten and molybdenum surfaces using a 


2 
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computer simulation model and to compare the results with the exper- 
imental data obtained by Winters. 

The incident ion will be Ar? in the energy~range from 0.5 keV to 
3.0 keV. 

The simulation includes: 

1) Bombardment of a clean (001) tungsten surface (used for 
comparison). 

2) Bombardment of a nitrogen reacted (001) tungsten surface. 

3) Bombardment of a nitrogen reacted (001) tungsten surface where 
the atomic mass of tungsten is replaced by the atomic mass of 
molybdenum. 

4) Bombardment of a nitrogen reacted (001) molybdenum surface, where 
the nitrogen is closer to the target surface plane than in the 
cases listed above. (The nitrogen is assumed to be in a 4-fold 


symmetry position for all reacted surfaces.) 


III. THE SIMULATION 


A. THE TARGET 
1. Properties of Materials 


a) Tungsten (W, Wolfram): atomic mss 183.85, lattice parameter 
3.165 A, crystal structure BCC, nearest-neighbor distance 


2.74 A, cohesive energy 8.9 eV, 


ie 





b) Nitrogen (N): atomic mass 14.0, binding energy to the tungsten 
substrate 6.5 eV, 

c) Molybdemun (Mo): atomic mss 94.94, lattice parameter 3.147 A, 
crystal structure BCC, nearest-neighbour distance 2.73 A, 
cohensive energy 6.82 eV. 

2. Location of the Adatoms 

The location of nitrogen atoms on a tungsten or molybdenum 
surface is not known. Clavenna [21] and Adams [22] suggest that the 
nitrogen is single-bounded in a 4-fold position above an atom in the 
second layer of the tungsten substrate. (The 4-fold position is also 
assumed by Winters.) Figures 2a,b show side and top view of the 
reacted surface where the radius of the nitrogen is assumed to be 0.59 
A {22]. (The smaller circles are the nitrogen atoms.) Criffiths 
et. al. [23], agree with this position at a fractional coverage (0Q¢) 
of 0.5 but mention, that at this coverage the nitrogen is unstable. 
They recommend study at the most stable coverage of O¢ = 0.4. This 
requires formation of islands of nitrogen atoms on the surface, (see 
Figs. 3a,b)- Harrison et. al. [24], propose that the nitrogen might be 
below the first layer of the substrate in an interstitial site. In 
this case some second layer atoms Would be pushed down into the 
substrate to set equal distances to all nearest-neighbors. For this 
Simulation the nitrogen is assumed to be in the 4-fold position for 
both substrate materials. For the nitrogen reacted tungsten surface 
the location above the substrate is chosen so that the distance to all 
five nearest-neighbors is equal. The distance above the surface is 


0.504 LU (where 1.0 LU = 1.58 A = 0.5 lattice parameters, see Fig. 4). 
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(Lattice units Were introduced in the program so that rescaling to a 
new target material can be done without major reprogramming and so that 
results from different substrates can be compared directly.) Runs also 
were made with the tungsten mass replaced by the molybdenum mass to 
examine the mss effect proposed by Winters [18], discussed in chapter 
2. On the nitrogen reacted molybdenum surface the nitrogen is placed 
closer to the surface layer (see Fig. 4). The distance above the 
surface is 0.243 LU (where 1.0 LU = 1.57 A = 0.5 lattice parametrs). 

3. The Target Size 

The inability of the lattice of target atoms to contain the 
entire cascade generated by the incident ion is called ' a failure in 
containment’. The size of the target determines the containment of the 
trajectory. If an increase in target size doesn't change the results 
of the computations, the trajectory is effectively contained for the 
purposes of the computations. (Unnecessarily large targets waste 
computer resources. ) 

The minimum target size for the ion energy of 0.5 keV was found 
to be 19x8x19, where the numbers represent the numbers of planes in the 
x,y and z directions. 40 nitrogen atoms (equivalent to a fractional 
converage of O¢ = 0.4) are placed as described above (see Fig. 5; 
solid lines mark the smaller target). The minimum target size for ion 
energies between 1.0 keV and 3.0 keV was found to be 23x8x23 with a 
fractional coverage of O¢ = 0.416, equivalent to 60 adatoms. The 
difference in coverage does not significantly affect the result. 
Winters [18], determined the coverage in his experiments by Auger 


Spectroscopy, and found that the Auger intensity is a nonlinear 
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function of coverage above 9 = 6 x 1014 atoms per cm*. Thus all 
his measurements were taken below this value. 
014 


The coverage chosen for the simulation is equivalent to 4.9 x l 


atoms per cm?. 


B. POTENTIAL FUNCTIONS 
1. In General 

A primary uncertainty in sputtering models and theory is the 
form of the interatomic potential functions. WN single potential 
function has yet been deduced which fits every sputtering experiment. 
The functions used for this simulation are Morse, and Moliere [25] 
joined with a cubic spline [24]. The Morse potential function is of 
the form: 

VG) —sOAamerD laze) — 2D. exp|—a(r—-r,) | 
Where D, is the well depth, a is a scale factor of the potential 
function, and ro is the equilibrium separation. For r > rp the 
Morse potential is attractive and for r < ro repulsive. 
The Moliere potential function is of the form: 
V(r) = (ZZ 9e*/r) (0.35 exp[-0.3r/a] + 0.55 exp[-1l.2r/a] 
+0.10 exp[-6.0r/a]) 

where 'a' is the Firsov screening length. The Moliere function is a 
pure repulsive potential function. 

The cubic spline is of the form: 

V(r) = Cy + C; xr + Cor? + Cor? 

Where the coefficients are determined by the program. The cubic spline 


is necessary to combine Morse and Moliere functions into a single 
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composite potential. For hard collisions, at close approach in the 
keV-range, the Morse potential bends towards the Y-axis and thus cannot 
be physically correct, because all potential functions must approach © 
as r > 0 like l/r. The cubic spline is inserted where the Morse and 
the Moliere functions overlap, thus the combined function is a Moliere 
to the left and a Morse to the right of the spline. In general the 
spline shouldn't be longer than 0.05 LU to get smooth potential and 
force curves. 
2. Crystal Potential Functions 

The composite potential functions are used for the bulk 
materials as Well as for the reacted surfaces. The Moliere function is 
a standard Moliere in each case. The Morse parameters have physical 
meaning: a, used as an adjustable Parameter, is chosen so that the 
slope of the Morse matches that of the Moliere function. For the 
substrate materials D, is chosen to generate individual atom energies 
in the center of the target equal to the cohesive energy, i.e., 8.9 eV 
for tungsten and 6.82 eV for molybdenum. For the reacted surfaces De 
has to match the binding energies between the nitrogen and the substate 
iee., 6.5 eV in each case. Because the potential function is truncated 
after one nearest-neighbor distance, rp is greater than the 
separation between atoms in the substrate materials. For the reacted 
surfaces, Tf) is the nearest-neighbor distance. 

The cubic spline is bounded by the values Ry and Rpe The 


table below shows the parameters for the crystals amd adatoms. 


ey, 
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W -W 
W -N 
Mo-Mo 
Mo-N 


Table l: 


a( Ant) 


1.20 
1.65 
1.50 
2. 20 


Parameters for Crystal Potential Function. 


aan) 


2-849 
2-380 
2-800 
1.950 


DaCev ) 


1.355 
1.300 
0.997 
ji. G2 


Ra(A) 


WSS 
1.14 
1.24 
0.91 


Rp (A) 


1.90 
1.14 
1.30 
0.97 


see Figs. 


O70 
i ao 
5.0 
9.10 


Figures 6, 7, 8, 9 show the Morse and Moliere functions and Figure 
10 shows the combined functions generated by the program. 
3. Gas Phase Substrate Atom Potential Functions 
The gas phase potential functions are pure Morse functions 


because no hard collisions are expected. The parameters are listed 


below: 

Table 2: Parameters for Gas Phase Substrate Atom Potential Functions. 
a(A71) re(A) DeCeV) see Figs. 

W —W 1.20 Dey a 8-900 11 

Mo-N ZeZ0 1. 96 5.100 eZ 


4. Nitrogen-Nitrogen Potential Functions 
The N-N potential function in the gas phase as well as in the 


crystal phase is a pure Morse functions with: 


a = 2.7 A7t 

Yo = 1.098 A 

Dalat 3 eV (crystal phase) 
Da = 9.834 eV (gas phase) 


see Figures ll, 12 for the gas phase potential functions. 


5. Argon-Tungsten/Molybdenum and Argon-Nitrogen Potential 


Functions 
Because it's not likely that argon ions react with nitrogen 
tungsten or molybdenum, these potential functions have the standard 


Moliere form, (see Fig. 13). 
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C. THE IMPACT AREA 

The impact area for each simulation is a rectangle containing a 
grid of 104 impact points, (see Fig. 14), where the length of the sides 
is 3.1654 = 2.0 LU for tungsten and 3.147 A = 2.0 LU for molybdenum. 


The origin of the rectangle is placed at the following locations: 


Table 3: Origns of the Impact Area 


target size locations 

23x8x23 N-reacted CLOeEO) (10,92) ie. 5 
19x8x19 Nereacted (08,08) (08,06) Bee 5 
23x8x23 clean Cie 1 1) Page 
19x8x19 clean (09,09) Fig. 5 


On the Nitrogen-reacted surface, two impact point Sets are neces- 


Sary to assure complete converage. 


IV. DISCUSSION AND RESULTS 


All output was analyzed by a separate analyzer program, AN83, which 
provides a numerical analysis, and by the program PLOTAN, providing a 


graphical analysis. The latter is discussed first. 


Ae GRAPHICAL ANALYSIS 
1. Spot Pattern 
The spot pattern plots show the polar and azimuthal angular 
distribution of the ejected atoms, (see Figs. 15, 16, 17). The 
distances of the crosses form the center is a measure of the polar 
angle of the ejected atoms, i-e., a cross far out corresponds to a 


large ejection angle measured from the surface normal. The spot 


ihe 
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pattern obtained for the substrate mterials (W, Mo, Mo-mass) show a 
four-fold symmetry, where 'Mo-mass' means, that in the program the 
tungsten mass was replaced by the molybdenum mss and all other 
parameters were unchanged. 

The spots at 45°, 135°, 255° and 315° point in the direction of 
the nearest~neighbors. Figure 15 shows the spot patter for molybdenum 
atoms ejected from the substrate after bombarding a nitrogen-reacted 
molybdenum surface by 2.0 keV argon ions. This figure is typical of 
all spot pattern obtained from substrate mterials. 

The spot pattern obtained from the adatoms show no symmetry. 
Figures 16, 17 are representative for ejected nitrogen atoms. Figure 
16, the case where the nitrogen is placed further above the surface, 
shows a heavier concentration of ejected atoms at low angles than 
Figure 1/7. 

2- Energy Distribution 

Figure 18 shows typical energy distributions for atoms 
sputtered from the substrate in the ion energy range from 0.5 keV to 
3.0 keV. The envelopes peak near 7.0 eV. 

The energy distribution of the sputtered nitrogen atoms depends 
on the location of the adatom. If the nitrogen is close to the surface 
layer, the distribution peaks at about 2.5 eV/atom, (see Fig. 19). If 
the nitrogen is further away, as in the case of the tungsten or 


molybdenum—mass substrate, the distribution is much broader, and the 


maximum is roughly 15 eV/atom, (see Fig. 20). 
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3. Energy Distribution >100 Fsec 

The energy distribution for particles ejected from the 
Substrate after 100 fsec (1.0 fsec = 1071> sec) is similar for all 
materials, but with the difference that the fraction for molybdenum or 
molybdenum—mass is less than for tungsten, (see Figs. 21, 22, 23). 
More low energy atoms are ejected after 100 fsec in the case of the 
heavier substrate material. This could indicate a mass effect. 

Energy distributions for the ejected nitrogen atoms are shown 
in Figures 24 and 25. If the nitrogen is in the lower position there 
were no atoms ejected with energy greater than 7.0 eV. This behavior, 
is seen throughout the whole energy range, (0.5 keV to 3.0 keV). 
Although the total number of atoms is low, (see Fig. 25), and therefore 
the statistics are poor, the difference between Fig. 24 and Fig. 25 my 
be indicating a distance effect. 

4. Energy Versus Time 

Figures 26, 2/7, 28 and 29 show that no high energy particles 
are ejected later than 100 fsec. The spectrum for substrate materials 
(Fig- 26) is wider than that for the ejected nitrogen atoms. 

Figure 28 shows that hardly any nitrogen atoms are ejected 
after 100 fsec. The difference between Fig. 27 and 29 and Fig. 28 on 
the other hand shows again, that there may be a distance effect. 

3» Ejection Time Distribution 

The ejection time distribution is essentially the same for all 
substrate materials. For high energies the distribution is lower and 
broader, (see Figs. 30 and31). The distribution peaks at about 80 


fsec.e 
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6. Atoms per Single Ion 
The atoms per Single ion diagrams are Similar throughout all 
runse The spectrum broadens with increasing energy, (see Figs. 38, 39, 
40 and 41). 
7. Yield per Single Impact Point 
The yield per single impact point diagrams are also similar for 
all runs, (see Figs. 42 and 43), where the open circles represent no 
yield. The diagrams are combined for the locations (10, 10) and (10, 
12), iee., the upper left corner of each rectangular box represents the 
location (10, 10). Comparison of the centers of the combined squares 
shows, that whenever the ions hit a substrate atom the yields are 
larger. The center of each right Square is occupied by a nitrogen 
atom. The nitrogen acts like a 'buffer'. The open circles are caused 
by channeling. Figures 42 and 43 show that with decreasing ion energy 
the areas where channeling occurs decrease. The reason for this is 
that for decreasing energies the cross sections increase because the 
impact energy appears in the denominator. 
8. Yield per Surface Atom 
Figures 44, 45 show the yield obtained from the surface layers 
for both target sizes. Almost no atoms are ejected from the edges of 
the crystal. This indicates that the crystals are large enough to 


contain the trajectories. 


Be NUMERICAL ANALYSIS 


The numerical analysis done by the analyzer program AN 83 provides 


the following nitrogen yields, (see also Fig. 46): 
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Table 4: Nitrogen Yields from Simulations. 


substrate: W(001) Mo-mass(001) Mo (001) 
energy yields yields yields 
in keV 

G5 2230 2.04 0. 60 
1.0 225 2-29 0.68 
Zoi O 1.80 269 0.51 
3.0 leo? 1.87 Ceol 


Experimental yields in the ion energy range from 0.5 keV to 3.0 keV 
are known only for polycrystalline tungsten and molybdenum surfaces 
[26]. Figure 47 shows the experimental metal yields obtained from 
polycrystalline tungsten and molybdenum as well as the yield from clean 


WC001) obtained from the simulation. The yields are listed below: 


Table 5: Yields from Clean Metal Surfaces. 


Experiment Simulation 
energy WC poly) Mo( poly) WCOO1) 
in keV 
0.5 0.58 0.57 2031 
sO 0. 88 0. 80 2.94 
2.0 oe) 1.10 2.40 
3.0 1.50 ie Pao 


Some data are known for 5.0 keV Ar’ bombardment of molybdenum 

[20] at normal incidence onto the (001) planes and onto polycrystalline 
molybdenum. The yield from polycrystalline Mo is found to be larger 
than from Mo(001). Thus Roosendaal [20] concludes that for a BCC 
crystal the same result is found as for FCC structures: the sputter- 
ing yield for less closed packed planes is less than the yield from 
polycrystalline targets, i-e., Y(Poly) > Y(001). The results of this 
Simulation show a yield of about twice the yield from polycrystalline 
tungsten. Winters [16] and Wehner [27] who bombarded W(001) surfaces 
in the low energy range with various ion types obtained the results 
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shown below: 


Table 6: Experimental Yields from Clean Metal Surface for Ion 
Energy of 0.5 keV. 


Neon Argon Krypton Xenon 
Winters: 0.25 a Omar ieee L 
Wehner ° On 2] ad OQ. 92 ve O2 


The yield is increasing with increasing ion mass and should be 
about 0.6 for 0.5 keV Ar* ions. The result of this simulation is 
higher than this value by a factor of about 4. This shows that the 
(Ar-W) potential (a standard Moliere) is too hard in this energy 


range. 
Winters [18] shows his results as cross section versus ion energy, 


(see Fig. 1). The nitrogen cross sections are: 


Table 7: Experimental Nitrogen Cross Sections. 


cross sections x 107}, (cm) 
energy substrate WC(001) Mo(001) 
in kev 
OF5 0.91 0.65 
1.0 0.95 O./2 
200 1.00 0.75 
cme 0.94 O73 


Winters assumes that the nitrogen sputtering rate is given by 
-Ry = dOy/dt = oNOnvt, 


Where Oy is the nitrogen atom concentration (atoms/cm), On 


+ the ion flux (ions/cm*- 


is the sputtering cross section (cm2), v 
Ss), and t the time in seconds. The product tv? is the ion dose 
(fluence) while the product 9,0n gives the sputtering yield in 


nitrogen atoms/ion. 
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Using this relation (Y = 0,0) the cross section obtained from 


the simulation are; 


Table 8: Simulation Nitrogen Cross Sections. 


cross sections x 19715, (em?) 
energy substrate: WCOO1) Mo-mass (001) Mo(001) 
im keV 
Oe5 Gad 4.42 lie 2ZO 
ae0 ans 4.63 lwd6 
200 3.64 ie. 42 iO? 
360 Qe] 7 3 uerS desi? 


Figures 48 and 49 show Winters cross Sections in comparison to the 
cross sections obtained from the simulation. in the case where the 
nitrogen is in the lower position (Mo) (Fig. 48), the simulation cross 
sections are in a better agreement with the experimental results. 3ut 
even in this case the Simulated cross sections are still too large. 
This suggests that the (Ar-N) potential (also a standard Moliere) is 
too hard. Winters' experiments showed that the sputtering yield of 
nitrogen tends to increases as the atomic weight of the substrate 
increases. He suggests that the difference in the cross section for 
the sputtering of nitrogen from W(001) and Mo(001) is dominated by the 
adsorbate~substrate mss difference since the other important 
parameters are similar. The simulation shows that the distance between 
the nitrogen and the substrate is more important than the mss 
dirtference. 

Winters furthermore States that the shape of the yield curves for 
the chemisorbed gas and the substrate material is different because the 
ratio o(gas)/o(substrate) decreases with increasing ion energy. 

Winters used o(substrate) cross sections obtained from polycrystalline 
Wand Mo assuming a surface coverage of 9 = 10715 atoms/cm*. The 


following ratios were calculated from the simulation: oN(W)/oW(clean), 
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o(Mo-mass)/oW(clean), oN(Mo)/oW(clean). Using oW(clean) in the 
latter case seems reasonable, because comparision of yields from 
polycrystal- line W and Mo as well as the yields obtained from the 
substrates are not too different, i.e., the yields from Mo(clean) 
shouldn't diviate much from W(clean). © for W(001) is 9.98 x 10+! 


atoms/cm. The ratios are shown below (see also Fig. 49). 


Table 9: Cross Section Ratios. 
Experiment Simulation 
oONCW) o N(Mo) o NCW) oN(Mo-mass) oN(Mo) 


energy oW(poly) oMo( poly) oW(clean) oW(clean) oW(clean) 


in keV 


0.5 ey) 1.14 2-06 1.78 0.52 
1.0 1.11 0. 90 1.54 ero 0.46 
2-0 0.77 0.68 1.52 1.84 0.43 
3.0 0.63 0.56 1.2/7 1.73 0. 47 


In comparison to Winters’ results, Figure 50 shows the ratios for 
nitrogen in the lower position and Figure 51 for nitrogen in the higher 
position. (Mo-mass not shown.) The values obtained from the simula- 
tion show the general trend observed by Winters for the W and Mo 
substrate, but indicate that the location of the high nitrogen is 
perferable for the tungsten substrate and low for the molybdenum 
substrate. Even With these choices, the difference between the 
nitrogen locations is too large. This is also seen by comparing 
oN(W)/o N(Mo) between Winters data and the data obtained from the 


simulation. The ratios are shown below, (see also Fig. 52): 


Table 10: Cross Section Ratios. 


energy Experiment Simulation 
in keV oNCW)/o N(Mo) oONCW)/oN(Mo-mass) oN(W)/o N(Mo) 
0.5 1.40 1.16 | 3.98 
Las) 1.36 0.98 3655 
2.0 33 O. 82 Deo 
3.0 1.29 0.73 Zaid 
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As mentioned in section A/7, with decreasing ion energy the cross 
sections increase and the areas where channeling occur decrease. This 


is also seen comparing the reflected argon yields, listed below: 


Table 11: Simulation Argon Yields. 


clean nitrogen reacted 
energy substrate  W(001) Mo-mass(001) W(00O1) Mo(001) 
in keV 
0.5 0.40 QO. 30 0.43 OQ, 24 
0 Os 7 Oeil Oran 2 0.15 
7.0 OeuZ Q.07 0.09 0.05 
Dre) 0) 0.07 QO. 07 0.05 0.03 


Due to the increase in cross sections, more lions are reflected at 

low energies. Harrison (19) says that about 85% of the ejected atoms 
(bulk material) from clean targets come from the first layer. This 
Simulation shows that at least 80% of the ejected atoms come from the 
first layer. The yield from deeper layers varies between 19% and 5%, 
where the yield from below the second layer is almost neglible. The 


yields from below the first one are shown below: 


Table 12: Yields from Second or Deeper Layers (Simulation). 


clean nitrogen reacted 
energy substrate WCOOL) Mo~mass(001) WC(001) #£Mo(001) 
in keV 
0.5 11% iy 16% 5% 
1.0 10% 14% er 8% 
2.0 16% 16% 19% 12% 
3.0 L7% 18% 17% 16% 


Some of the ejected atoms, which have similar energy, velocity and 
position, combine above the surface, forming multimers. Most of the 
multimers observed are dimers: (W-W, (W-N), Mo-Mo), (Mo-N). No (N=N) 


dimers were found. The numbers of dimers observed are listed below: 
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Table 13: Dimers from Clean Tungsten (Simultion). 


energy W-W 
in keV 

0.5 40 
1.0 53 
2.0 18 
3.0 D2 


Table 14: Dimers from Nitrogen Reacted Surfaces (Simulation). 


energy Wo-N Mou- N W- W Mo - Mo 
in keV 
0.5 50 9 18 lis 
1.0 65 2 23 Sill 
2.0 Z| 15 24 15 
3.0 14 8 2 10 


The table below shows the number of linked atoms (taken from all 


multimers). 


Table 15: Linked Atoms (Simulations). 


energy substrate W Mo-mass Mo 

in keV (WW) CWN) (MoMo) (MoN) (MoMo) (MoN) 
0.5 35 rae 56 85 33 18 
1.0 41 82 74 a7 54 ao 
2.0 i 66 oF 64 27 18 
3.0 42 Bi 42 65 29 i! 


In the case of the molybdenum substrate the number of Mo-N links 
exceeds the number of Mo-Mo links throughout the whole energy range. 
In the case of the tungsten substrate the number of linked W-W is 


higher for 3 keV and 2 keV. 
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Ve SUMMARY AND RECOMMENDATIONS 


The results of the simulation show that the ion-substrate 
potential functions are chosen too hard. Further work should be done 
to reduce the yields for the substrate materials. The graphical 
analysis shows that there may be mass as well as distance effects, but 
the results of the simulation clearly show that the adsorbate-substrate 
distance is more important than the adsorbate-substrate mass 
difference, when the nitrogen sits in the higher position on the 
substrate with the higher atomic mass. After adjusting the substrate 
yields, futher work should be done to find the right distance between 


the adatoms and the substrate. 
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